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ABSTRACT 

We present a detailed analysis of the XMM-Newton RGS high resolution X-ray spectra 
of the Seyfert 2 galaxy, Mrk 573. This analysis is complemented by the study of the 
Chandra image, and its comparison to optical (HST) and radio ( VLA) data. The soft 
X-ray emission is mainly due to gas photoioniscd by the central AGN, as indicated 
by the detection of radiative recombination continua from O vii and O viii, as well 
as by the prominence of the Ovii forbidden line. This result is confirmed by the 
best fit obtained with a self-consistent cloudy photoionisation model. However, a 
collisionally excited component is also required, in order to reproduce the Fexvii lines, 
accounting for about 1/3 of the total luminosity in the 15-26 A band. Once adopted the 
same model in the Chandra ACIS data, another photoioniscd component, with higher 
ionisation parameter, is needed to take into account emission from higher Z metals. The 
broadband ACIS spectrum also confirms the Compton-thick nature of the source. The 
imaging analysis shows the close morphological correspondence between the soft X-ray 
and the [O ill] emission. The radio emission appears much more compact, although 
clearly aligned with the narrow line region. The collisional phase of the soft X-ray 
emission may be due to starburst, requiring a star formation rate of ~ 5 — 9 Mq yr^^, 
but there is no clear evidence of this kind of activity from other wavelengths. On the 
other hand, it may be related to the radio ejecta, responsible for the heating of the 
plasma interacting with the outfiow, but the estimated pressure of the hot gas is much 
larger than the pressure of the radio jets, assuming equipartition and under reasonable 
physical parameters. 
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1 INTRODUCTION 

Important progress has been made in the last few years to 
unveil the origin of the soft X-ray emission in obscured Ac- 
tive Galactic Nuclei (AGN). The first breakthrough was rep- 
resented by high resolution spectra made available thanks to 
the gratings aboard Chandra and yiMM- Newton. The 'soft 
excess' observed in CCD spectra was found to be due to 
the blending of bright emission lines, mainly from He- and 
H-like transitions of light metals and L transitions of Fe, 
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with low or no continuum, in most Seyfert 2 galaxies (see 
e.g. Sako et al. 2000; Sambruna et al. 2001; Kinkhabwala 
et al. 2002; Guainazzi & Bianchi 2007). Spectral diagnostic 
tools agree that the observed lines should be produced in a 
gas photoionised by the AGN, with little contribution from 
any collisionally ionised plasma. A second breakthrough was 
made possible thanks to the unrivaled spatial resolution of 
Chandra. The soft X-ray emission of Seyfert 2 galaxies ap- 
pears to be morphologically correlated with that of the Nar- 
row Line Region (NLR), as mapped by the [Olll] A5007 
HST images (e.g. Young et al. 2001; Iwasawa et al. 2003; 
Bianchi et al. 2006). Since the NLR is also believed to be 
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a gas photoionised by the AGN, it was shown that a very 
simple model where the soft X-ray emission and the NLR 
emission are produced in the same material is possible (e.g. 
Bianchi et al. 2006). 

However, this scenario is clearly oversimplified. In par- 
ticular, it is not clear whether different components of the 
same medium are spatially separated, possibly radially dis- 
tributed, or co-exists at each radius, as a result of stratifi- 
cation. Moreover, the role of radio ejecta, which arc often 
found to be strongly correlated with the morphology of the 
NLR, is still unclear. An exciting possibility to shed some 
more light on this issue is provided by spatially resolved 
X-ray spectroscopy, which is, unfortunately, limited by the 
compact structures of this class of objects (at most some 
arcsec). Indeed, such a study was performed with high res- 
olution spectroscopy on only one source, NGC 1068, being 
very bright and extended (Brinkman et al. 2002). The re- 
sults are in agreement with the expectations from a cone of 
plasma, irradiated by the central AGN. On the other hand, 
a similar analysis, but with CCD resolution, was performed 
on NGC 7582, showing that there are regions with a further 
source of ionisation, or lower density (Bianchi et al. 2007). 

Mrk 573 (a.k.a. UGC 1214, z=0.0172) is the third [O ill] 
brightest source in the Schmitt et al. (2003) sample of 
nearby Seyfert galaxies observed by HST, only fainter than 
NGC 1068 and Mrk 3. Moreover, the NLR extension of 
Mrk 573 is by far the largest of the sample, reaching a total 
extent of 3 kpc, with a projected extension almost 9 arc- 
sec wide. A triple radio source is associated to the gala:xy, 
composed by a central core and two spots (Ulvestad & Wil- 
son 1984). A detailed analysis of the NLR of this source 
was performed, among others, by Ferruit et al. (1999) and 
Schlesinger et al. (2009). They conclude that photoionisa- 
tion by the central AGN is likely the dominant process, al- 
though the interaction with the radio jets must be taken 
into account in the overall scenario, introducing kinematic 
disturbances and shaping the NLR morphology. 

Despite the brightness and the interesting features of 
Mrk 573, the source has never been studied in detail in X- 
rays. After the detection by Einstein, Mrk 573 was observed 
by ROSAT, which clearly detected an emission extended 
over about 10 arcsec, in agreement with the NLR dimen- 
sions. The source was then only observed by XMM-A^euiion 
with a short exposure time of about 10 ks. It confirmed to 
be rather bright in the 0.5-2 keV band (~ 3 x 10"" cgs) 
and the detection of a very strong iron line, with an EW 
larger than 1 keV, revealed its nature as a Compton-thick 
source (Guainazzi et al. 2005). But the most important piece 
of information comes from the RGS high resolution soft X- 
ray spectrum, which clearly appears dominated by strong 
emission lines. The predominance of K lines, the ratio of the 
components of the O vii triplet and the detection of strong, 
narrow, radiative recombination continua (RRC) features all 
contribute to an interpretation in terms of photoionised gas 
(Guainazzi & Bianchi 2007). 

In this paper, we re-analysod in detail the XMM-Newton 
RGS spectrum of Mrk 573, adopting a self-consistent pho- 
toionisation model. Moreover, we present for the first time 
the imaging and spectral analysis of a Chandra observation. 



2 OBSERVATIONS AND DATA REDUCTION 

In the following, errors correspond to the 90 per cent con- 
fidence level for one interesting parameter (Ax^ = 2.71), 
where not otherwise stated. The adopted cosmological pa- 
rameters are Ho = 70 km s"'^ Mpc~^, Qa = 0.73 and 
Q,n = 0.27 (i.e. the default ones in xspec 12.5.1: Arnaud 
1996). At the distance of Mrk 573, 1 arcsec corresponds to 
360 pc. In all the fits, the Galactic column density along the 
line of sight to Mrk 573 is included (Nh = 2.96 x 10^° cm"^: 
Dickey & Lockman 1990). 

2.1 X-rays: Chandra and XMM-Newton 

Mrk 573 was observed by Chandra on 2006-11-18 for a to- 
tal exposure time of 40 ks (obsid 7745), with the Advanced 
CCD Imaging Spectrometer (ACIS: Garmire et al. 2003). 
Data were reduced with the Chandra Interactive Analysis 
of Observations (CIAO: Fruscione et al. 2006) 4.1 and the 
Chandra Calibration Data Base (CALDB) 4.1.2 software, 
adopting standard procedures. Images were corrected for 
known aspect offsets, reaching a nominal astrometric ac- 
curacy of 0.6 arcsec (at the 90 per cent confidence level). 
Spectra were extracted from three different regions: a circu- 
lar region of 16 arcsec of radius (the default one analysed in 
Sect. 3.2); the same region, but excluding the inner 1 arc- 
sec; the rmclear region, i.e. the inner 1 arcsec. The imaging 
analysis was performed on event files without the pixel ran- 
domization and treated with the SER procedure (Li et al. 
2003), in order to improve the positional accuracy. This al- 
lowed us to use a pixel size of 0.246 arcsec. 

Mrk 573 was also observed by XMM-A''ewton on 
2004-01-15 for a total exposure time of ~ 12 ks (obsid 
0200430701), with the EPIC CCD cameras, the pn (Striider 
et al. 2001) and two MOS (Turner et al. 2001), operated 
in Prime Full Window and Medium Filter. Those data were 
already presented by Guainazzi et al. (2005) and Guainazzi 
& Bianchi (2007). In this paper, we present a new and more 
detailed analysis of the RGS spectra, which were extracted 
with standard procedures, using SAS 8.0.1 (last described 
in Gabriel et al. 2004) and the most updated calibration 
files available at the time the data reduction was performed 
(October 2009). Background spectra were generated using 
blank field event lists, accumulated from different positions 
on the sky vault along the mission. 

As shown in Fig. 1, in the field there are several 
bright sources in the soft X-ray band. However, the ex- 
traction region (corresponding to the default 90 per cent 
of the point spread function in the cross-dispersion di- 
rection) includes only one contaminating source, named 
[TUM93] J014401.5-h022106 (Turner et al. 1993). We ex- 
tracted its pn spectrum, and calculated a 0.5-0.8 keV flux of 
1.6 X 10"" erg cm"^ s"^, i.e. a factor of 8 less than Mrk 573. 
It is unlikely to contaminate the line emission spectrum of 
Mrk 573, since any line, if present, would be displaced by 
the order of 1 A (without considering the effect of its un- 
known redshift), given its distance of around 1 arcmin from 
the nucleus. Such lines would therefore be at unidentified 
wavelengths, but none of them are detected (see Sect. 3.1). 

We also extracted the pn spectrum of Mrk 573, already 
analysed by Guainazzi et al. (2005). However, since the qual- 
ity of the data of this short observation is lower than that of 
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Figure 1. XMM-Newton EPIC pn image (0.3-2 keV). The RGS 
extraction region is shown in green. 

the Chandra ACTS spectrum, we will not discuss it further 
in this paper. 

2.2 Radio: VLA 

We downloaded from the VLA archive the observation of 
Mrk 573 performed on 1985-03-03 at 6 cm (4860 MHz) in A 
array configuration. Data reduction was performed with the 
NRAO software package AIPS, following standard proce- 
dures, and very little flagging of bad data points. The phases 
of the calibrator 0146-1-056 were interpolated and applied to 
the object, and 3C 48 was used as the primary flux cali- 
brator. The final map was produced using the AIPS task 
IMAGR with a beam size of 0'.'5x0'.'4, cleaning depths of 
several thousand iterations and a CLEAN gain factor of 0.1. 

2.3 Optical and IR: HST 

We retrieved the HST observations we use in this paper from 
the Multimission Archive at STScI (MAST). The images 
were processed through the standard on-the-fly reprocess- 
ing system. Mrk 573 was observed with HST with WFPC2 
and the FR533N ramp fllter as part of the GO program 
6332 on 12/11/1995. The target was located in the WF2 
camera to image the [Olll]A5007 emission line. The total 
exposure time for the two CR-SPLITted exposures is 600s. 
We combine the two images using the crrej task on IRAF 
to allow for cosmic-ray rejection. The image is clearly dom- 
inated by the emission line gas, therefore for the purposes 
of this work, a careful subtraction the continuum emission 
is not needed. The object was also observed in the near- 
IR with HST/NICMOS, as part of G07867. Images were 
taken with both the FllOW and F160W filters, which are 
similar to the J and H-bands, respectively, using the NICl 
camera. The exposure time is 1023s for both images. The 
data were originally presented by Martini & Pogge (1999). 
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Table 1. Mrk 573: detected emission lines and RRC in the XMM- 
Newton RGS spectrum. 
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in eV. Theoretical energies are from CHIANTI (Dere et al. 1997, 
2009). The labeUing for Fexvii Unes follows that of Brown et al. 
(1998). 

The FllOW filter includes the Pa/3 emission line, while the 
F160W filter is relatively free of strong emission lines, at 
the redshift of Mrk 573. The image at longer wavelengths is 
clearly dominated by the unresolved AGN emission. In order 
to assess the possible presence of structures in the central 
regions of the galaxy (e.g. starburst regions, dust), we sub- 
tract the stellar emission of the host galaxy using a model 
obtained by fitting ellipses to the galaxy isophotes. The re- 
sult shows that no extra nuclear structures are present. This 
image was therefore used as a reference for all the other im- 
ages, by aligning the IR nucleus to the soft X-ray and [O in] 
brightest pixels. 



3 SPECTRAL ANALYSIS 

3.1 The soft X-ray RGS spectrum 

Following the procedure described in Guainazzi & Bianchi 
(2007), the simultaneous fits on the spectra of the two RGS 
cameras were performed on ~ 100-channel-wide segments, 
adopting the Cash statistics (Cash 1976). All the fits include 
a power law component (F fixed to 1^) and as many emis- 
sion lines as required (at the 90 per cent confidence level). A 
total of ten emission lines and two RRC were detected and 
identified securely (see Table 1). The widths of all the fea- 
tures were consistent with being unresolved. The RRC are 
then properly modelled with the redge model, which fully 
takes into account the recombination edge profile, giving a 
direct estimate of the electron temperature of the gas. 

The resulting temperature (< 15 eV) is much less than 
would be expected from the same species, if the gas were 
in coUisional equilibrium. This is a strong piece of evidence 
in favour of a dominant photoionised phase in the emitting 

1 Given the very limited band of these fits, the modellisation of 
the continuum is insensitive to the power law photon index, which 
can be frozen at any value. 
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Figure 2. Mrk 573: centroid energy vs. flux contour plots for the detected emission lines in the XMM- Newton RGS spectrum (see Table 
1). The contours refers to AC=2.30, 4.61 and 9.21, i.e. confidence levels of 68, 90 and 99 per cent for two interesting parameters. Upper 
limits for the closest lines are also plotted (AC=1.00, 2.71 and 4.00). 
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region. However, the Ovil RRC (0.7393 keV) may be sig- 
nificantly contaminated by the 3F component of Fcxvil L 
emission (0.7390 keV). This is suggested by the detection of 
the blend of the 3G and M2 components of the same species, 
at ~ 0.7266 keV. Our simulations with the APEC model in 
XSPEC showed that the 3F/(3G-|-M2) ratio never exceeds 
~ 0.44, for a wide range of temperatures. On the other hand, 
the observed ratio between the Une detected at 0.7389 keV 
and that at 0.7266 keV is larger, being 1.2'tQ'l. This means 
that a large part of the observed flux must be due to Ovil 
RRC. The O viii RRC may also be contaminated by Fe xviii 
L lines (see Table 2), but no other emission lines from that 
ion arc detected in the spectrum. 

In principle, diagnostics on the Ovil triplet may also 
be a good indicator of the ionisation mechanism of the gas. 
Indeed, the predominance of the forbidden component is a 
sign of photoionisation, and the significant detection of the 
resonant transition a hint that pure recombination is not the 
only mechanism to produce the emission lines, but photoex- 
citation has also an important role. This, in turns, suggests 
that the gas column density should be not too large, in or- 
der not to suppress this process. However, given the large 
uncertainties on the line fluxes, it is impossible to exclude 
other solutions only on these grounds, like a contribution to 
the resonant line from a coUisionally excited plasma. 

Another hint for photoionisation comes from the Ovil 
Ka forbidden line to Oviii Ka ratio, which is larger than 1 
{2.5't\'f,)- Together with a large total luminosity of the Ka 
oxygen lines of 6 x 10*" erg s~^, this puts Mrk 573 in the 
photoionisation-dominated locus of the empirical diagnostic 
plot presented in Guainazzi et al. (2009). 

Since the phenomenological analysis of the RGS spec- 
trum favours an origin in a photoioniscd gas, we tried a more 
physical approach, performing a self- consistent fit on the 
whole spectrum (limited by the S/N to the 15 — 26 A band). 
We produced a grid model for xspec using cloudy 08.00 
(last described by Ferland ct al. 1998). The main ingredients 
are: plane parallel geometry, with the flux of photons strik- 
ing the illuminated face of the cloud given in terms of ion- 
isation parameter U (Osterbrock & Ferland 2006); incident 
continuum modelled as in Korista et al. (1997); constant 
electron density Uo — W'' cm^"^; elemental abundances as 
in Table 9 of cloudy documentation^; grid parameters are 
log?7 = [-0.25 : 2.00], step 0.25, and logA^H = [19.0 : 23.5], 
step 0.1. Only the refiected spectrum, arising from the il- 
luminated face of the cloud, is taken into account in our 
model. We also produced tables with different densities 
(ue = 10^ — lO'' cm~^): all the fits presented in this pa- 
per resulted insensitive to this parameter, as expected since 
we are always treating density regimes where line ratios of 
He- like triplets are insensitive to density (Porquet & Dubau 
2000). 

The fit with a single photoionised phase is presented 
in the upper panel of Fig. 3. With logU = 0.6 ± 0.3 and 
log Nh — 20.6 ±0.7, a good fit is obtained, with most of the 
lines (ictcctod in our phenomenological analysis reasonably 
modelled, including, notably, the RRCs. However, there are 



^ Hazy 1 version 08, p. 47: http://viewvc.nublado.org/index. 
cgi/tags/release/c08 . 00/docs/hazy 1_08 . pdf ?revision= 
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clear residuals at the Fexvil 3G+M2 wavelength: these lines 
are completely missing from the photoionised model^. This 
is not surprising, given the very low ionic fraction of Fe XVII 
(~ 4 X 10~*) at this ionisation parameter. 

We therefore tried to add another photoionised phase, 
with a larger U. However, any attempt to reproduce the ob- 
served Fexvii lines largely overpredicts the Oviii Ka and 
K/3 lines. This forced us to abandon this approach and to 
try an alternative scenario, where the iron L lines are mainly 
produced in a coUisional gas. Indeed, the addiction of a col- 
lisional phase (kT = 0.30lo;Qg keV) perfectly takes into ac- 
count the Fe xvii 3G-I-M2 lines, together with the 3C com- 
ponent at 0.826 keV, without affecting too much the rest of 
the spectrum. The improvement of the Cash statistics for 
the addition of the coUisional gas component is AC = 18, 
with two less degrees of freedom. No further improvement is 
achieved allowing the elemental abundances to vary. The pa- 
rameters of the photoionised phase change to log U = O.llJJ ^ 
and logNn = 20.515 • The residuals of this hybrid (pho- 
toionised-|-collisional) gas are presented in the lower panel 
of Fig. 3 and are rather good. The only exception is rep- 
resented by some positive residuals for the Ovii K7, but 
they are only significant in the RGSl spectrum. Indeed, a 
separate analysis of the two RGS spectra leads to a line de- 
tection only in the RGSl (2.5 ± 1.4 x lO"'^ ph cm"^ s"^ at 
0.6995 ± 0.0013 keV, formally inconsistent with O VII K7), 
while in the RGS2 only an upper limit is measured if the 
line energy is fixed at the theoretical one (< 0.8 X 10-^ ph 
cm~^ s~^). We therefore conclude that the RGSl detection 
is insecure, and therefore cannot draw any conclusions based 
on that only*. 

The total flux in the 0.5-0.8 keV band is 1.3 x 10~" 
erg cm^^ s^^. The contribution of the coUisional phase is 
4.6 X 10"" erg cm"^ s"\ i.e. roughly 1/3 of the total flux 
in this band. 

3.2 The Chandra broadband spectrum 

The broadband (0.4-8 keV) Chandra X-ray spectrum of 
Mrk 573 appears dominated by a strong emission in the 
soft band, plus strong neutral iron Ka emission, typical sig- 
natures of an highly obscured AGN. The high energy part of 
the spectrum is well fitted by a pure neutral reflection com- 
ponent plus the iron line. The properties of the reflection 
component are unconstrained, due to the low statistics, so 
wc fixed the photon index to 1.7 (as in typical Scyfert galax- 
ies, see e.g. Bianchi et al. 2009) and the cosine of the incli- 
nation angle to 0.45. The largest normalization compatible 
with the data leads to an iron EW of ~ 1.8 keV, perfectly in 
agreement with the expectations for a Compton-thick AGN 
(Matt et al. 1996). Since the normalization of the reflection 
component is basically unconstrained in the following fits, 
where the soft X-ray component may contribute also to the 
high energy part of the spectrum (see below), we decided to 

The c;loudy line database was rather inaccurate for the Fe xvn 
transition wavelengths and relative atomic parameters. We there- 
fore modified it according to the CHIANTI database. 
* We note here that the line wavelength does not cor- 
respond to any known defective pixel in neither of the 
RGS cameras (http://xnim.esac.esa.int/external/xnim_user_ 
support/documentation/uhb/node59.html#3177). 
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Figure 3. XMM-Newton RGS spectra (15 — 26 A) of Mrk 573, rebinned for displaying purposes only. Top: Best fit with a single 
photoionised phase, whose parameters are reported above. Bottom: Best fit with a photoionised and a collisional phase, whose parameters 
are reported above. The brightest emission lines for the adopted models are labelled. 



fix it. The 2-10 keV flux and the iron hne flux are in per- 
fect agreement with the values measured with yiMM-Newton 
(Guainazzi et al. 2005). 

To flt the soft X-ray band, at first we adopted the 
same model which successfully reproduces the RGS spec- 
trum. This model is a good representation of the ACTS spec- 
trum, when limited to the same band, but clearly fails to fit 
the remaining part of the data, because it cannot reproduce 
the observed emission from higher Z metals (see left panel 
of Fig. 4). We added another photoionised phase to the fit, 
with a larger ionisation parameter, and allowed also the pa- 
rameters of the other phases to vary. The resulting flt is very 



good (x^ = 74/79 d.o.f., see right panel Fig 4 and Table 2). 
The best fit model includes two further emission lines, which 
can be readily identified with neutral Si Ka and S Ka at 
1.740 keV and 2.308 keV, respectively (House 1969), likely 
arising from the same Compton-thick material responsible 
for the production of the features which dominate the high 
energy part of the spectrum, that is the Compton refiection 
component and the neutral iron Ka line. 

The best fit parameters for the two photoionised phases 
are logU = O-Stoi and logNn = 21.7toj, and logU = 
1.8ll|^:i| and logNn = 21.8toi- The temperature of the 
collisional gas is kT = 0.56 ± 0.12 keV. 
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Figure 4. Left; Mrk 573 Chandra ACIS-S broadband spectrum, with tiie best fit for tlie RGS data in tlie soft X-ray part, and Ax^ 
residuals. Right: The same, but with the best fit. 



Table 2. Best fit parameters for the Chandra spectrum of 
Mrk 573. Fluxes are in units of 10"^"^ erg cm~^ s~^, line fluxes 
in units of 10~® ph cm~^ s~^, energies and kT in keV. See text 
for details. 



logUi 
log Nhi 

logU2 
log Nh2 

kT 

EsiKci 

FsiKa 

EsKa 

FsKa 

EpeKa 

FpeKa 

xVdof 

^0.5-2 
^2-10 



0.3: 

21.72 
1.8lj 

21.8"* 



-■-0.6 
7-1-0.6 
-0.8 
,-H0.15 
-0.12 
i+0.7 
-0.8 

0.56 ±0.12 
1.740* 
1.0 ±0.7 
2.308* 
1.2 ±0.8 
6.35 ± 0.03 
5.5 ±1.7 
74/79 

2.7 ±0.3 

2.8 ±0.5 



fixed 



If the three-phase model is re-applied to the RGS spec- 
tra, the fit is visually comparable to the one obtained above 
directly on these data, with only two phases. Therefore, 
although only 2 phases (one photoionised, the other col- 
lisional) is required by the RGS data, the Chandra ACIS 
data reveals the presence of a further photoionised com- 
ponent, still allowed by the XMM-Newton high-resolution 
spectrum. 

The total 0.5-2 keV flux is 2.7 x 10"^^ erg cm"^ s"\ 

in agreement with the one measured with X.MM-Newton 
(Guainazzi et al. 2005). The contribution from the coUisional 
phase is ~ 20%, corresponding to an unabsorbed luminosity 
of 4.0 X 10^° erg s~^ in the same band. It is interesting to 
note that, of the remaining photoionised phases (which con- 
tribute roughly equally to the the soft X-ray flux), ~ 15% 
is constituted by a continuum component (which includes 
RRCs and the Thomson-scattered power law), while the re- 
maining flux is in emission lines. 



Table 3. Flux (mjy) and luminosity (erg s ^ Hz ^) at 6 cm for 
the three radio components detected in Mrk 573. 

Flux logL 

core 0.95 ± 0.03 27.72 

NW 2.22 ±0.03 28.09 

SE 4.54 ± 0.05 28.40 



4 IMAGING ANALYSIS 

The soft X-ray (0.2-2 keV) emission of Mrk 573 is clearly ex- 
tended and closely resembles the morphology of the NLR, as 
mapped by the [Olll] HST image (see left panel of Fig. 5). 
This is often observed in Seyfert 2 galaxies, and clearly sug- 
gests a common physical origin for the two emissions (e.g. 
Bianchi ct al. 2006). Given the lower angular resolution of 
Chandra, it is hard to say how good the correspondence 
between the X-rays and the NLR is at the smallest scale. 
However, there is no evidence that significant deviations are 
present in these data. 

On the other hand, the right panel of Fig 5 shows the 
VLA radio emission superimposed on the Chandra data. As 
already reported by Ulvestad & Wilson (1984) and Falcke 
et al. (1998), Mrk 573 has a unresolved radio core of a di- 
ameter size of ~0.32 kpc (1 arcsec) and two-sided jets: the 
northwest and the southeast bubbles have a similar size of 
-0.48x0.32 kpc^ (l'.'5xl'.'0). Table 3 reports fiux and lu- 
minosities for these three radio components. It is clear that, 
although the inclination of the radio jets is aligned with that 
of the soft X-ray emission, the overall extension of the radio 
emission is far more compact. 

We also tried to do some spatially resolved spectroscopy, 
to look for variations of the properties of the soft X-ray emit- 
ting gas along the distance from the nucleus. Following the 
analysis performed on the RGS high resolution spectra, we 
chose two narrow energy bands, which we know are domi- 
nated by emission from Ovii (0.5-0.6 keV) and Oviii (0.6- 
0.7 keV)^. As shown in Fig. 6, there is an hint that the O Vii 

^ By means of simulations with the observed RGS fluxes of the 
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Figure 5. Left: Chandra soft X-ray (0.2-2 keV) contours superimposed on the HST Olll image. The contours refers to 0.001, 0.1, 0.2, 
0.4 and 0.95 levels with respect to the brightest pixel. North is up, east to the left. Right: VLA radio (6 cm) contours superimposed on 
the Chandra soft X-ray (0.2-2 keV) image. 
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Figure 6. Mrk 573: Chandra images in the 0.5-0.6 kcV (Ovil) and 0.6-0.7 keV (Ovill) band. The images were smoothed with a 3x3 
FWHM Gaussian filter and normalised to the brightest pixels. They are re-scaled logarithmically in the same way (see colourbar at the 
top). The black cross is the brightest pixel of the hard X-ray image. North is up, east to the left. 



emission lines and the spectral response of the Chandra observa- 
tion, we estimated that the 0.5-0.6 keV band is contaminated by 
less than 10% from Ovill photons, while the 0.6-0.7 keV band by 
around 25% from O VII photons. 



emission is more extended than the O vill one, but the qual- 
ity of the data does not allow us to quantify this difference. 

We also extracted two spectra, one only of the inner 1 
arcsec region, the other excluding it, while keeping all the re- 
maining X-ray emitting region. While practically all the flux 
above 2 keV is concentrated in the nuclear region, approx- 
imately 15% of the 0.5-2 keV flux is produced farther than 
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1 arcsec from the nucleus. This percentage is significantly 
lower than what is generally found in Seyfert 2s observed by 
Chandra (see e.g. Bianchi et al. 2006). The best fit param- 
eters for the three phases of the soft X-ray emission in the 
nuclear region are consistent with those found for the outer 
region. There is an hint of an higher relative flux between 
the collisional and the photoionisation gas in the outer re- 
gion with respect to the nucleus, but it is not statistically 
significant. 



5 DISCUSSION 

5.1 Fluorescence from Compton-thick material 

The high energy spectrum of Mrk 573 is typical of a 
Compton-thick Seyfert 2 galaxy. While the Compton reflec- 
tion component and the strong iron Ka line axe ubiquitous 
signatures of reflection from a Compton-thick material, the 
presence of fluorescent Ka lines from neutral silicon and 
sulphur are less common in spectra of obscured AGN (but 
see e.g. Sambruna ct al. 2001; Bianchi et al. 2005; Tilak 
et al. 2008). This is due to the fact that the fluorescence 
yield is a strong function of the atomic number, and is thus 
much lower for low-Z metals with respect to iron. Using the 
Basko (1978) formulae, valid in the case of a semi-infinite 
plane-parallel slab isotropically illuminated, we calculated 
the expected ratios between the fluxes of the silicon, sul- 
phur and iron Ka fluorescent lines. We adopted fluorescence 
yields of 0.042, 0.077 and 0.304, respectively (Kaastra & 
Mewe 1993), and the Anders & Grevesse (1989) abundances. 
For all the lines, we considered both the unscattered and 
the once-scattered photons, since with the quality of our 
data we cannot disentangle the two components even in the 
iron fine. The calculated ratios are Si/Fe=0.010 — 0.028, 
S/Fe=0.014 - 0.032, Si/S=0.70 - 0.88, where the ranges 
take into account the dependencies on the inclination angle 
and the incident power law index. These values, once trans- 
formed in EWs with respect to the incident continuum, arc 
very similar to the ones calculated by Matt et al. (1997). 

The observed ratios are Si/Fe=0.18±0.14, S/Fe=0.22± 
0.16, Si/S=0.8 ± 0.8. The fluxes of the observed sificon and 
sulphur fluorescent lines arc significantly higher than pre- 
dicted, although the uncertainties are quite large. It is dif- 
ficult to recover the observed fiuxes with elemental over- 
abundances, which would be very high, and apparently not 
required by the ionised plasma dominating the soft X-ray 
emission. It is likely that the observed fluxes of these two 
lines are rather unreliable, because of the low statistics and 
resolution of the spectra. Indeed, the neutral Si Ka is proba- 
bly contaminated by the Mgxii fine, at 1.745 keV, which 
could be underpredicted by our photoionisation model. We 
note that the fluxes of neutral Si and S Ka are also large in 
other sources where they are tentatively detected (e.g. Sam- 
bruna et al. 2001; Bianchi et al. 2005; Tilak et al. 2008). 

5.2 The soft X-ray emitting region 

The soft X-ray emission of Mrk 573 appears dominated by 
a photoionised gas, morphologically coincident with the op- 
tical NLR, as commonly found in Seyfert 2 galaxies (e.g. 
Bianchi et al. 2006; Guainazzi & Bianchi 2007). From the 



Chandra images, we can see that the two, symmetric, emit- 
ting regions have an inner radius < 90j5C (pixel size of 
0.246 arcsec), and an outer radius ro — 1 kpc. The pro- 
jected opening angle of the cones is roughly 50 degrees, in 
agreement with the 45 degrees measured by Wilson & Tsvc- 
tanov (1994). Since the intrinsic luminosity of this source 
can only be indirectly estimated, because the nucleus is ob- 
scured by a Compton-thick material, it is difficult to derive 
a better guess for the inner radius of the emitting region, 
from the ionisation parameters of the fits with photoionisa- 
tion models. 

On the other hand, from the luminosity of the strongest 

line, the O vil forbidden line (which is not significantly con- 
taminated by the collisional phase), wc can calculate its 
emission measure (EM), and estimate the density of the ma- 
terial producing it. The luminosity (ph s"'^) of a recombina- 
tion line of charge state i is: 

Li= Az.U+irja{T)nenHdV (1) 
Jv 

where Az is the abundance of the element, /j+i is the charge 
state fraction of the recombining ion, a{T) is the radiative 
recombination coefficient for the recombining ion at temper- 
ature T, 77 is the fraction of recombinations leading to the 
relevant transition, and Ue and hh are the electron and hy- 
drogen density in the volume V of the emitting region (e.g. 
Liedahl 1999). The EM is defined as: 

EM^ [ n,n„ dV = , (2) 
Jv Azfi+iva{T) 

where the last equivalence assumes that all the relevant 
parameters are constant in the gas (or their average val- 
ues arc used) . Adopting the best- fit parameters of the RGS 
fit, CLOUDY gives the following values: Az = 4.9 x 10"", 
fi+i = 0.13, a{T) = 1.3 X 10"" cm^ s"^ (for an average 
temperature of 4.7 x 10'' K), and rj = 0.48. From the ob- 
served luminosity of the O vil forbidden line, we get an EM 
of 9 x 10^^ cm~'^. The volume of the emitting region may be 
approximated (r, <^ ro) to F ~ g[l — cos(7r/4)]7rro/6, where 
g is the filling factor of the emitting gas within the spheri- 
cal bi-cone with opening angle of 45 degrees. By assuming 
Ue ^ 1.2nH, we derive: 

gril ~ 25 cm~^ (3) 

Therefore, the filling factor g is likely to be much smaller 
than 1. If we take the average densities measured in the op- 
tical NLR (of the order of 10^ - 10^ cm"^: Capetti et al. 
1996; Ferruit et al. 1999; Schlesinger et al. 2009), we de- 
rive g ~ 10~^ — 10~®, i.e. the emitting gas fills a very small 
fraction of the whole volume of the bi-conical region. Such 
a small filling factor is also required in order to recover a 
column density of the order of 3 x 10^° - 6 x 10^' cm"^, 
measured from the RGS and ACIS spectra, along the ob- 
served ~ 1 kpc. This result is in agreement with the very 
low mean density (~ 1 cm~^), derived from the EM of H/S, 
if a unity filling factor is assumed (Schlesinger et al. 2009). 

5.3 The collisional phase and the radio emission 

The detailed spectral analysis of the RGS spectra of Mrk 573 
pinpointed the presence of a gas phase in collisional equilib- 
rium, which contributes by ~ 20 — 30% to the overall soft 
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X-ray luminosity. The observed thermal emission could be 
in the hot gas surrounding a starburst region. If all the 0.5-2 
keV luminosity emitted by the collisional phase is associated 
to star formation, we can estimate a star forming rate (SFR) 
of ~ 9 M0 yr"^ (Ranalli ct al. 2003). A SFR of the same 
order (~ 5 Mq yr~^) can be estimated from the total radio 
luminosity of the two spots, if related completely to the star 
formation, assuming the relation found by Bell (2003), and 
an average spectral index of a = —0.75 (Falcke et al. 1998). 
A lower SFR (~ 2 Mq yr^^) is derived from the far in- 
frared (FIR) emission (where the starburst emission largely 
dominates over the AGN, e.g. Fritz et al. 2006) adopting 
the relation found by Kcnnicutt (1998), and the FIR lu- 
minosity calculated from the IRAS infrared fluxes (Moshir 
et al. 1990) and the method proposed by Helou et al. (1988). 
Moreover, as already reported in Sect. 2.3, there is no clear 
evidence of star-forming regions in the inner few kpc in the 
HST ncar-IR images. 

On the other hand, the collisional phase may be directly 
connected to the observed radio emission, both arising as 
free-free emission of a hot gas. The X-rays to radio ratio for 
thermal bremsstrahlung emission of a plasma at temperar- 
ture T is: 

(4) 

g{vr,T)e "Ft 

where g {v, T) are the Gaunt factors appropriate for the X- 
rays {vx) and radio (vr) frequency, at that temperature (e.g. 
Longair 1992). Adopting a temperature of 5 x 10® K, of 
the order of the one measured for the collisional plasma in 
the RGS and ACIS spectra, the Gaunt factor at 5 GHz is 
around 11, while at 1 keV is approximately 1.5 (Longair 
1992). Therefore, the luminosity ratio between the X-rays 
and the radio free- free emission should be ~ 0.01. For a radio 
luminosity of 10^*^ erg s^"*^ Hz^"*^ (representative of each of 
the regions detected in the VLA data, sec Table 3), we would 
expect an X-ray luminosity of 10'^*' erg s^^ Hz^^ at 1 keV, 
i.e. 10"'^ '' erg s~^ keV~^, several orders of magnitude larger 
than what is observed. 

Note that this docs not exclude that all the radio emis- 
sion that we observe in Mrk 573 is due to free-free emission 
of a hot plasma. The nuclear core, for example, is likely to 
be absorbed by a substantial column density of neutral gas, 
like the X-ray nucleus, so that the X-ray emission possi- 
bly associated with the radio emission would be completely 
suppressed. On the other hand, the radio spots, extended on 
larger scales, should not be affected by such a large obscu- 
ration, and cannot be due to thermal bremsstrahlung emis- 
sion, without overproducing the observed X-ray emission by 
a large factor. In any case, it is clear that the collisional 
plasma we detect in the X-rays cannot be produced by a 
thermal plasma which gives rise also to the radio emission. 

The radio ejccta may still be responsible for the heat- 
ing of the X-ray-emitting plasma, though not emitting in 
X-rays themselves. In this scenario, the radio luminosity of 
the two jets is mainly due to non-thermal synchrotron, while 
the steep power law index (Falcke et al. 1998) makes their 
X-ray luminosity negligible. Indeed, Capetti et al. (1996) as- 
signed a crucial role to the radio ejecta, because they com- 
press the line-emitting gas, enhancing the emission where 
this interaction occurs. This interaction could also be the 



heating source of the gas at ~ 5 x 10*^ K that we observe 
through its signatures of plasma in collisional equilibrium. 
This scenario could be tested by making a comparison be- 
tween the pressure of the hot X-ray emitting gas and the 
minimum pressure of the radio jets. The latter can be esti- 
mated by assuming equipartition between the particles and 
the magnetic field. We model each jet lobe as a sphere of 
radius 0.2 kpc. We assume the electron energy distribution 
extends from a Lorentz factor 7min = 2 up to 7max = 10^, 
with an electron energy index p = 2.4 (the results have only 
a weak dependence on the choice of 7; see HardcEistle et al. 
2004). Under the simple assumption that the jet is in the 
plane of the sky, and that relativistic beaming is unimpor- 
tant, the minimum pressure of each lobe is 3x 10^^^ barye. 
On the other hand, the pressure of the hot gas can be es- 
timated from the density directly derived from the APEC 
normalisation of our Chandra best fit. Assuming the same 
conical emitting volume as in the previous Section, but with 
an outer radius of 500 pc, similar to the one observed for 
the radio emission, we get a pressure of ~ 1.4 x 10~^ barye 
(a very similar value, ~ 1.1 x 10~^ baxye, is derived if we 
use the APEC normalisation of our RGS best fit). 

The pressure of the hot gas is therefore about 50 times 
larger than the minimal internal pressure in the radio jets. 
For the shocks to be the source of heating of the gas, the 
two pressures should be comparable. In order for the jet not 
to be suppressed (unless we are not seeing it at some special 
time), there must be some additional pressure. Indeed, it 
was suggested that departures from equipartition may char- 
acterise the weak jets and lobes observed in Seyfert galax:ies 
and FR I-type radio galaxies (see e.g. Evans et al. 2006, 
for the case of NGC 2110), differently from the powerful 
jets observed in FR Il-typc radio galaxies. However, other 
studies of the NLR in Mrk 573 have shown that the the opti- 
cal spectrum is well reproduced by nuclear photoionisation, 
and there is no evidence of shocks produced by the radio 
outflows (Ferruit et al. 1999; Schlcsinger et al. 2009). The 
latter do not likely have any strong influence on the NLR 
medium, apart from some kinematic effects due to their ex- 
pansion into the gas, eventually influencing the formation of 
the observed arcs (Schlesinger et al. 2009). 



6 CONCLUSIONS 

We presented a self-consistent analysis of the XMM-Newton 
RGS spectra of the Seyfert 2 galaxy, Mrk 573. Several pieces 
of evidence suggest that the dominant ionisation process of 
the soft X-ray emitting gas is photoionisation: the clear de- 
tection of RRC from O vil and O vill, and the prominence 
of the O VII forbidden line. A photoionisation model fully 
takes into account all the brightest emission lines, but a 
collisional phase is also required, in order to reproduce the 
Fexvil lines. This component accounts for about 1/3 of the 
total luminosity in the 15-26 A band. 

The broadband Chandra ACIS spectrum confirms the 
Compton-thick nature of the source, dominated by a Comp- 
ton refiection component and a strong neutral iron Ka line. 
The soft X-ray emission needs a further photoionisation 
component with respect to the RGS spectrum, with a larger 
ionisation parameter, in order to reproduce emission from 
higher Z metals. 
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The Chandra soft X-ray image closely follow the NLR 
morphology mapped by the [Olll] emission. On the other 
hand, the radio emission is far more compact, although 
clearly aligned with the NLR. It could be directly related 
to the coUisional phase found in the X-ray spectra, in a 
plasma heated by the interaction with the radio ejecta, but 
the estimated pressure of the hot gas is much larger than 
the pressure of the radio jets, assuming equipartition and 
under reasonable physical parameters. Alternatively, the gas 
in coUisional equilibrium may originate in a starburst region, 
requiring a star formation rate of ~ 5 — 9 Mq yr~^, but there 
is no clear evidence of this kind of activity from other wave- 
lengths. Deeper X-ray observations arc needed in order to 
confirm the presence of a gas in coUisional equilibrium in 
Mrk 573, and understand its nature. 
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